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ABSTRACT 

The timing properties of the 4.45 s pulsar in the Be X-ray binary system GRO J 1750-27 are 
examined using hard X-ray data from INTEGRAL and Swift during a Type II outburst observed 
during 2008. The orbital parameters of the system are measured and agree well with those 
found during the last known outburst of the system in 1995. Correcting the effects of the 
doppler shifting of the period, due to the orbital motion of the pulsar, leads to the detection 
of an intrinsic spin-up that is well described by a simple model including P and P terms of 
-7.5 X 10"'" s s"' and 1 x 10 s s"^ respectively. The model is then used to compare the 
time-resolved variation o f the X -ray flux and intrinsic spin-up against the accretion torque 
model of Gh osh & Lamb! (1 19791) : this finds that GRO J1750-27 is likely located 12-22 kpc 
distant and that the surface magnetic field of the neutron star is ~ 2 x 10'^ G. The shape of the 
pulse and the pulsed fraction shows different behaviour above and below 20 keV indicating 
that the observed pulsations are the convolution of many complex components. 

Key words: pulsars: individual: GRO J1750-27 - stars: neutron - X-rays: binaries - X-rays: 
bursts - ephemerides. 



1 INTRODUCTION 

The number of known high mass X-ray binary systems (HMXB) 
has expanded in recent years with the advent of hard X-ray satel- 
, lite telescopes with large fields of view enabling the serendipitous 
' detection of these sources whilst in outburst. In particular the num- 
, ber of HMXB systems comprising a neutron star (NS) with a Be 
star companion, which in most cases are only seen in brief out- 
bursts following ~ years of quiescence, has increased to the point 
where the majority of HMXB known now are such systems (here- 
after BeXRB). Be stars are young and hot, with spectral Type B 
(or O) characterised by Balmer emission lines thought to be asso- 
ciated with circumstellar material shed due to the rapid rotation 
of the star, which may be near to its disruption speed (see e.g. 
IPorter & Riviniusll2003h . The NS, in an elliptical orbit around the 
star, can accrete matter from an equatorial "decretion disk" that 
may be formed from the disrupted material, or from a strong stellar 
wind; this potentially leads to outbursts in X-rays seen at different 
orbital phases. 

The types of ou tburst can be separated into Type I and Type 
IIl lStellaetai]|l986l) : Type I outbursts are seen to occur around the 
periastron of the binary orbit (although outbursts are not necessar- 
ily seen at every periastron) and are thought to be due to an increase 
in accretion onto the NS's surface when it is closest to the Be star 
where density of matter available is higher; Type II outbursts are 
much brighter than Type I, do not seem to be correlated with or- 
bital phase and in many cases last for several orbital periods. Pulsed 
X-rays are often seen from BeXRBs since they contain strongly 



magnetised NS; to date there are no known BeXRBs containing a 
confirmed black hole. Typical pulsar periods, f , for BeXRBs are in 
the range 1-100 sees a nd show a cor relation with the binary periods, 
P„b, of 10-100 days Corbetll 19861) . During Type II bursts the pul- 
sation rate of the pulsar is seen to increase {P reduces, or frequency 
V = l/P increases) suggesting that significant angular momentum 
transfer from the star to the NS is happening; this phenomena is 
known as spin-up. 



The observation of Type II outbursts are particularly interest- 
ing for many reasons, especially those where a pulsar can be ob- 
served: with luminosities of 1-2 xlO""* ergs s"' they are a factor 
of ~ 10 brighter than Type I outbursts and should allow good sig- 
nal to noise even for systems located far away from Earth; the be- 
haviour of the system over at least one full orbital period can be 
investigated, allowing the orbital charactersitics to be inferred from 
changes in P\ the matter accretion rate, m, can be studied by relating 
the spin-up rate, P, of the pulsar to the transfer of angular momen- 
tum from the Be star; a period of spin-down of the pulsar towards 
the end of a large outburst with large spin-up may be indicative that 
the pulsar has undergone a torque reversal caused by an interaction 
of the maximally spinning pulsar magnetosphere with the inner ac- 
cretion disk - this would allow an estimate of the magnetic moment 
of the NS material to be examined; changes in the orbital phase 
averaged spectra may allow the Be star environment to be probed. 
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Figure 1. GRO J 1750-27 X-ray light curves from the Swift jBKi transient 
monitor webpage I http://swift.gsfc.nasa.gov/docs/swift/results/transients i 
and RXTE/ ASM (http://xte.mit.edu/ASM.lc.htmll, shown by vertical 
crosses and triangles respectively. Dark squares are overlaid on the BAT 
data to mark where simultaneous observations were made by INTEGRAL. 
The dashed hne shows the 15-50 keV flux from the X-ray burster GX 3+1, 
which is located 18' from GRO J1750-27 (average flux from 5 daily mea- 
surements is shown for clarity). 



1.1 GRO J1750-27 

The transient X-ray pulsar GRO J1750-27 (also known as 
AX J1749. 1-2639) was discovered towards the Galactic C entre, by 
the CGRO/BACrSE satellite experiment (Scott et al.lll997'), follow- 
ing an outburst of 20-70 keV photons over the period MJD 49915- 
49978 (1995 July 17 - September 18). The BATSE detectors were 
large wide-angle uncollimated scintillator detectors, hence the out- 
burst was detected only through a pulsed flux with period 4.45 s. 
Any un-pulsed component was not measurable, in the absence of 
any imaging capability, and formed part of the background sig- 
nal (along with all other sources in the FoV). The pulsed flux was 
equivalent to ~ 30 mCrab at its peak. The pulsation was observed 
to spin-up during the outburst at a rate correlated with the pulsed 
intensity, reaching 38 pHz s"' at the peak of the outburst. A mod- 
ulation of the p eriod was f o und to correspond to a 29.8 day orbital 
period and led lScott et J] ( Il997h to identify GRO J 1750-27 as a 
BeXRB un dergoing a gia n t (Type 11) outbur st following the clas- 
sification o f ICorbetl ( [T98^ . IScott et al.l ( Il997l) made an estimate of 
the distance to the source of at least 18 kpc from Earth. Due to the 
non-imaging nature of the BATSE detectors this measurement is 
based on an assumed pulsed fraction of 30% for GRO J 1750-27. 

A second outburst was reported by the Swift/BAT tran- 
sient monitor during an observation starting at 11:26:40 UT on 
29th Jan 2008 (MJD 54494). The source was first detected at 
0.031 ± 0.005 cts/cm^sec (140 mCrab) in the 15-5 keV en- 
ergy range of the BAT detector ( Krimm et al. 2008). The out- 
burst is shown in Fig. [T] and consists of a sharp rise to the 
peak flux in ~ 30 days and a slower decay over the following 
~ 120 days. The XRT soft X-ray telescope, also onboard Swift, has 
given the currently most accurate position for GRO J 1750-27 as 
RA(2000) = 1 7''49°'13U,Dec. (2000) = -26°38'37'/3 with an un- 
certainty of 5" tearthelmv et al4i2008.fl Note that the location of 



GRO J 1750-27 is ~ 18' from the bright X-ray burster GX 3+1, 
which requires that observations must be made with instruments 
with good angular resolution. No optical counterpart brighter than 
R ^ n mag has been reported, although the assumed distance to 
GRO J1750-27 of 18 kpc and its location of (l,b) = (2° 4, +0^5) 
place it some distance beyond the Galactic Centre, where optical 
observations are very difficult due to absorption. 

This outburst was also observed by the ES A X-ray/y-ray satel- 
lite INTEGRAL during Ga lactic Bulge monitoring observations 
jKuulkers et al.l2008l , l2007h . The monitoring programme consisted 
of ~ 12 ksec observations of the Galactic Bulge repeated approx- 
imately once per orbital revolution of INTEGRAL (the orbital pe- 
riod of INTEGRAL is ~ 72 hours) and continued throughout the 
brightest part of the outburst. Figure [2] shows images taken with 
the INTEGRAL instruments during the peak of the outburst. The 
JEM-X instrument has finer angular resolution than ISGRl and 
hence GRO J1750-27 and GX 3+1 are clearly resolved in the 
low energy image. Nevertheless the two sources are also resolved 
with ISGRl and, due to the relatively softer spectrum of GX 3+1, 
GRO J1750-27 is many times brighter at higher energies. 

The all sky monitor (ASM) data from the RXTE satellite 
were also inspected to give further coverage of the outburst in 
the 2-10 keV range (see Fig. [T). Individual ASM observations 
of GRO J1750-27 have low signal to noise ratios, since the low 
energy flux was approximately a third of that at higher enrgies 
and because the measurements were based on typically only a few 
hundred seconds of exposure; the observations were also sparcely 
spaced, with several days between each. This meant that it was 
necessary to rebin the ASM light curve in order to reduce the un- 
certainty on the flux measurements to an acceptable level. Finally, 
the ASM data are affected by contamination from nearby sources, 
therefore this data was not considered further. 

Following the 2008 outburst, the data from INTEGRAL and 
Swift were obtained, allowing this fascinating source to be revis- 
ited with the advantage of hard X-ray data from simultaneous ob- 
servations with imaging capability and the benefit of long quasi- 
simultaneous exposures. 



2 DATA 

2.1 INTEGRAL 

The INTEGRAL satellite is a X-ray/y-ray obs ervatory that was 
launched by ESA on 17 October 2002 (Wi nkler et alj|2003h . It 
was designed with several instruments to offer fine spectroscopy 
and imaging in the 3 keV-10 MeV range. Two instruments in par- 
ticular are of interes t for this work: abov e 18 keV, IBIS/ISGRI 
dUbertini et al] l2003l : iLebrun et al] |2003) provides hard X-ray 
imaging in a large (30°) field of view with angular resolution ~ 12', 
point source location accuracy ~ 3', sensitivity ~ 18 mCrab for 
a 5<T source se en in an 1 .8 ksec observation in the 20-60 keV 
range (e.g. Kui ilkers et aD[2007); the JEM-X soft X-ray monitor 
fcund et alhOO^ consists of two identical coded-mask gas detec- 
tors sensitive to X-ray photons in the range 3-25 keV within a 
5° field of view and an angular resolution of 30". Only the JEM-Xl 
detector was operational during the observations described here. 



' A position of G RO J1750- 2 7, me asured with the ROSAT X-ray tele- 
scope, is given in IScott et aj] h997l) as RA = 17''49"'12»7, Dec. = 



-26°38'36'' (72000) with 4" accuracy, and is cited as Dennerl and 
Trumper, A&A, 1997, (in press); this paper, or any later one with the same 
authors, was not found in a literature search. 
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Figure 2. Images of the GRO J 1750-27 region taken with INTEGRAL dur- 
ing revolution 654 (2008 Feb 20; MJD 54516) in the 3-10 keV range 
with JEM-Xl (top) and 18^0 keV with ISGRl (bottom). The grid shows 
Galactic coordinates while the grey-scale represents statistical significance. 
White crosses mark the respective source positions. In the JEM-X image 
GRO J 1750-27 is detected with a peak significance of 36a" (correspond- 
ing flux 133 mCrab); the peak of GX 3+1 is 60o- (235 mCrab). In the IS- 
GRI image GRO J 1750-27 is detected with a peak significance of 117cr 
(194 mCrab); the peak of GX 3+1 is 13a- (21 mCrab) 



INTEGRAL makes observations of the Galactic Bulge, ap- 
proxi mately every 3 days d uring the monitoring of the Galactic 
Bulge jKuulkers et alj|2007h and is ideal for following up ~ month 
long bursts. The coded mask instruments JEM-X and ISGRI have 
been used to measure the pulsations from GRO J1750-27 in the 
range 3-70 keV. The data were reduced with the standard Offline 
Standard Analysis software versi on 7 (OSA?]! provided b y the Inte- 
gral Science Data Centre (ISDC: ICo"urvoisier et al.l2003h . 

Both ISGRI and JEM-X have good timing charactersistics, de- 
pending on the brightness of the source in question. Investigation 
of the 4.45 s pulse from GRO J1750-27 is easily within the capa- 
bilities of both INTEGRAL instruments. 



2.2 Swift 

The Swift/BKr instrument on the Swift satellite jGehrels et al.l 
|2004| ; iBarthelmv et alj|2005h is a wide angle coded mask system 
in a Low Earth Orbit, designed primarily to detect Gamma Ray 
Bursts (GRB). As such it has a very large field of view, approx- 
imately 2 sr, in order to maximise its chance of detecting a ran- 
domly occurring GRB. This means that during the course of its 
normal operations about 15% of the sky is observed and so most 
areas of the sky have some exposure during each 90 minute orbit. 
In survey mode, images are produced in histogrammes of approxi- 
mately 5 minutes length, allowing a light curve of the outburst to be 
produced with excellent time coverage as seen in Fig.[T] Of course, 
the survey image histogrammes are of no use for timing investiga- 
tions on time-scales less than the histogramme duration. However, 
1 second count rate data from the whole detector plane is available, 
which allows strong periodic signals of more than a few seconds 
length to be observed. In this respect, the observations from BAT 
are comparable to those made with BATSE in 1995; strong pulsed 
signals should be detectable, although no measurement of absolute 
fluxes will be possible from the rate data. 



3 INTEGRAL SPECTRUM AT THE PEAK OF THE 
OUTBURST 

Exisiting X-ray instruments are unable to use reflection or re- 
fraction in optical materials to focus photons with energies more 
than a few keV (note, however, that several new missions are 
planned that aim to achieve this in the future with various methods; 
Simbol-X, NuStar, NeXT, GRI). The instruments onboard INTE- 
GRAL rely on coded masks to cast a shadow on a position sensi- 
tive detector, which can then be used to reconstruct the direction 
of the incident light. Some contamination of the flux measurement 
for a particular object is therefore unavoidable in observations of 
fields containing many sources. Spectral analyses of objects in the 
Galactic Centre will therefore be challenging and especially so for 
GRO J1750-27 given its proximity to, and the relative brightness 
of, GX3+1 (seeFig.|2j. 

In particular, the 0SA7 tools for spectral analysis with JEM-X 
are optimised for a single on-axis source and for fields contain- 
ing many sources it is recommended to make spectra from fluxes 
taken from images made in discrete energy ranged This method 
is thought to be more immune to contamination of the spectra than 
the standard method, although it is not practical to make images in 
bands less then ~ 1-2 keV wide. A comparison of spectra made 
by both methods, for the peak of the outburst during revolution 
654, showed that they were consistent above 5 keV; at lower ener- 
gies there was a clear overestimate in the flux using the standard 
0SA7 spectral extraction method. 

For higher energy spectral analysis with ISGRI the detrimental 
effects of other sources are somewhat reduced since, in revolution 
654, GRO J1750-27 is brighter than GX 3+1 above 18 keV by 
about a factor of ten and is also the brightest source in the whole 
field of view by a factor of two. 

The broadband spectrum during the peak of the outburst was 
therefore calculated in the 5-70 keV range and fit with a cut-off 
power law; 



^ 0SA7 is available from the ISDC at |http://isdc.unige.ch] 



^ See Known Issues in the JEM-X analysis cookbook at 
|http://isdc.unige.ch/?Soft+download| 
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A(E) = k[^ 



IkeV/ 



(1) 



In order to investigate the energy spectrum of GRO J1750-27, 
the data from JEM-Xl and ISGRI were fitted simultaneously with 
a phenomenological model, namely a cut-off power-law function. 
The spectrum can be seen in Fig. [3] Systematics of 3% were added 
to JEM-X and 2% added to ISGRI and an inter-instrument cali- 
bration factor, defined as a factor applied to the ISGRI spectral 
response with lEM-X held constant, was also included as a free 
parameter. This gave a chi-squared value of Idof = 28.6/22 
and returned values of F = -0.15 ± 0.3 and = 6.0^^ J keV re- 
spectively (90% confidence) and an instrument calibration factor of 
1.49. Assuming that most of the energy output of the source is in 
the range 0.1-100 keV gives a total flux of 6.5x 10"' ergs cm"- s"', 
although the uncertainty on this is ±30% due to the variation on the 
measurement of F alone. 

A physical model, ascribing the spectral behaviour of the 
source to the inverse Compton scattering of soft photons by a ther- 
mal plas ma of elec t rons, was also fitted to the data. The COMPTT 
model o f iTitarchukI ( [1991 returned a reduced chi-squared value of 
1.24 and the following values: the plasma temperature 4.6±0. 1 keV; 
the plasma optical depth 6.4^q '; the the soft photon temperature 
was fixed at 0. 1 keV as it lies far outside of the range covered by 
our spectral data. The inter-instrument calibration factor was 1.40. 
This model follows essentially the same form as the cut-off power- 
law, as can be seen in Fig.|3] and suggests that the X-ray spectrum 
is due to the scattering of soft X-ray photons in an optically thick 
plasma. 

In both cases, the calibration factor returned was not equal 
to unity. This shows that there are some systematic effects in the 
spectrum that are not fully understood when dealing with moder- 
ately bright sources in crowded fields. Since the fields of view of 
ISGRI and JEM-X are not the same size, the backgrounds of both 
instruments will be dependent on different objects. The objects con- 
tained in both fields of view will also not have the same effects on 
each instrument, due to the large range of the X-ray spectra of the 
objects in the Galactic Bulge; each source may have a different rel- 
ative brightness to GRO 11750-27 in the energy bands of the two 
instruments, and will therefore be of changing importance to both 
background estimations. 

The detailed analysis of the GRO 11750-27 spectrum is 
clearly complicated even when the source is at its brightest; under- 
standing the instruments' spectral response during the observations, 
and hence the spectral behaviour of GRO 11750-27 throughout the 
outburst and its variation with pulsar phase, will be the subject of 
future work. In order to investigate the spectral behaviour, with- 
out model dependent spectral fits, simple spectral hardness ratios 
have been calculated from the total count rate fluxes measured by 
INTEGRAL during the outburst: 



{F, - Fj) 
(F, + Fj) 



(2) 



where F/ is the 20-50 keV ISGRI count rate and and Fy is the 5- 
20 keV JEM-X count rate. Figure|4]shows the time evolution of the 
broad band spectral hardness in the 5-20 and 20-50 keV ranges for 
the total (pulsed+unpulsed) data. There is a gradual softening of 
the total flux from GRO J 1750- 27 during the outburst, with some 
indications of more complex variability. 




Energy (keV) 

Figure 3. The INTEGRAL spectrum of GRO J1750-27 during the peak of 
the outburst in revolution 654. Measurements with JEM-Xl cover the 3- 
20 keV range (vertical crosses), whilst ISGl^ data is taken from 18-50 keV 
(diagonal crosses). The lines show models which have been simultaneously 
fitted to both data sets as described in section[3] comptonisation (continuous 
line); cut-off power-law (dashed-dotted line). 
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Figure 4. INTEGRAL measurements of the hai'dness ratio, ^fi+fji . 
F, is the 20-50 keV ISGRI count rate and and Fj is the 5-20 keV JEM- 
X count rate. The dashed vertical lines show the times of periastron of the 
GRO J 1750-27 system. 



4 HARD X-RAY TIMING ANALYSIS 

4.1 ISGRI pulsation detection at 18-40 keV 

The 0SA7 tools iijight and barycent were used to create 18^0 keV 
light curves, with 0.5 s binning, and to correct the time stamps of 
the results for the orbital motion of the INTEGRAL satellite, for 
each set of Galactic Bulge monitoring observations. 

The ISGRI light curves analysed here are made from several 
fixed pointings separated by slews and, as is common in high- 
energy astrophysics observational analysis, may contain irregularly 
distributed time data af fected by telemetry gaps e tc. Therefore the 
Lomb-Scare le method jLombl 1976l ; IScarglell982l) . which is robust 
against these effects, is used to search for periodic signals. The er- 
ror on the angular frequency of any detected pulsation is given by 
the formulations of lHome &Baliunasl ( [l98^ : 
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where the amplitude A is calculated thus, 



N 



(3) 
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Table 1. 18^0 keV pulsations detected from INTEGRALflSGRl observa- 
tions ofGRO J1750-27. 



(Tn is the standard deviation of the noise in the lightcurve, 
estimated to be the same as the standard deviation of the lightcurve, 
0"j, for a weak pulsed signal; zo is the power of the pulsation from 
the periodogramme; N is the number of data points; T is the total 
length of the data set. 

In each of the 17 Galactic Bulge monitoring observations the 
Lomb-Scargle periodograms showed a significant peak at a fre- 
quency of v~0.22 Hz. This corresponds to a period of ~4.45 s 
which is in agree ment with the observed spin period identified by 
IScott et all ( Il997b . Table [T] lists the detected pulse period and the 
associated error for each observation. The ISGRI data clearly show 
periodic signals which can be measured with very high accuracies, 
typically ~ few parts in 10^. 



4.2 ISGRI measurements of pulsar spin-up 

The GRO J1750-27 pulse period, P, as a function of observation 
time is shown in Fig.[5]and clearly indicates that the observed spin 
period evolves with time. While the oscillatory motion may be ex- 
plained by a doppler shifting of the spin-period by the binary orbit 
of the system there is also an additional component resulting in an 
intrinsic spin-up of the pulsar 

The observed data was fit with a simple spin-up model for the 
pulsar given by: 



P(t) 



P(to) + P{t-to)- 



Pit-toY 



(5) 



This model is convolved with a standar d orbital mod el assuming 
the orbital parameters measured by Sco tt et aP ( 1 19971) . excluding 
the orbital period. This combined model when fit to the data re- 
sulted in a chi-square value of ;t'-/dof = 12.78/13. An orbital period 
of 29.806 ± 0.001 days was found which is compatible with that 
previously observed; there is also a clear measurement of P and 
P. The details of the fit parameters and the orbital parameters of 
IScottetalJ ( fl997l) are given in Table[2] The resultant best fit model 
is presented by the solid line in Figure[5] 

Refitting the data without assuming any of the parameters of 
IScott et alj ( fl997l) yields an acceptable fit. However, whilst the or- 
bital parameters are compatible with the previous observation they 
have larger associated errors; therefore the prefered fit is based on 
the original parameters. Removing the P parameter results in an 
unacceptable fit, ;>'^/dof = 50.04/14. An ftest on these fit results 
shows that the addition of P is highly significant, with a chance im- 
provement probablility of only 3 x 10"'. Removing both the P and 
P terms returns an even poorer fit by many orders of magnitude: 
X^/dof ~ 6000/15. Hence the intrinsic spin-up model of equation |5] 
is determined with a very high level of confidence. 

4.3 Pulsation measurements with JEM-X 

Following similar methods to the ISGRI analysis, data from the soft 
X-ray monitor onboard INTEGRAL, JEM-X 1, were reduced with 
OSA7 for a broad 5-20 keV energy band. The analysis of the JEM- 
X data returned values of P consistent with the ISGRI results for 
all contemporaneous observations, albeit at lower statistical signif- 
icance. In several of the earlier observations significant power was 
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Table 2. The orbital parameters for GRO J 1750- 27 derived from lScott etal] 
i 19971) and from the analysis of ISGRI 20-40 keV data. 
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Figure 5. GRO J1750-27 orbital solution based on ISGRI measurements 
(crosses). Note in most cases the error bars on the period are much smaller 
than the size of the crosses. 
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also seen at ~ 2.22 s, corresponding to half of the ISGRI period. In- 
spection of the phase-folded JEM-X light curves showed that below 
20 keV the pulse amplitude was weaker and in some cases showed 
a double peaked profile. This is investigated in section l431 

4.4 Searching the Swi/if/BAT data for a periodic signal 

The ~ 2 sr field of view of BAT means that many observations 
are made of the source location; if suitable power spectra can be 
produced then the spin history will be well covered in time. Since 
the survey mode data, from which images can be produced, are 
integrated over ~ 5 mins they are not usable for timing studies 
of GRO J 1750- 27. The only BAT data with suitable time res- 
olution is the detector plane count-rate, in which the total flux 
on the whole detector plane is read every second. Since no posi- 
tional information is available from these data, any pulsed signal 
from GRO J1750-27 is contaminated by the signals from all other 
sources in the very large field of view. This can be effectively con- 
sidered as an increase in the background noise which will appear 
non-statistical. No time dependent background correction is avail- 
able for the count-rate data in the current Swift standard analysis 
and consequently it has not been possible to extract reliable pulsed 
fluxes. The BAT does, however, have the advantage of many obser- 
vations of GRO J1750-27 throughout the outburst. 

Data were obtained for all pointings between MJD 54480 
and 54640 (2008 January 15 and June 23) where the position of 
GRO J1750-27 was within 30° of the pointing direction. This was 
to obtain a sub-set of data not too badly affected by noise towards 
the edge of the detector plane and so that the source was detected 
with a reasonable sensitivity. The data were broken up into 50 ks 
segments, which were defined consecutively from the start of the 
observation, and were each searched for a periodic signal. Due to 
the low Earth orbit of Swift and its pointing strategy the data are 
generally arranged in a few short sections of ~ 1000 seconds dura- 
tion within the time bin. 

To improve the chances of detecting the pulsar spin period a 
Bayesian statistical approach is applied to the data set. The gener- 
alised Lomb-Scargle periodogram method of Bretthorst (2001) is 
used with the period range of interest restricted to the narrow range 
4.445^.460 s i.e. tightly bounding the periodicity detected by IS- 
GRI. 

Figure |6] shows the resulting periodograms displayed as a 
function of time and coloured for the probability density at each pe- 
riod. The highest probability density corresponds to periods which 
are consistent with the ISGRI measured pulses and the fitted period 
model; this is shown by the central dashed line in Figure |6] Also 
evident are side lobes of lower probability which follow the same 
time evolution; these are offset from the main track by integer units 
of a constant, k = ±3.4 ms. If this effect is caused by the beating of 
the pulsation with another periodic signal is given by the following: 



P P + k 

Inserting P = 4.450 s gives P(j„, ~ 97 minutes, which corresponds 
approximately to the orbital period of the Swift satellite. The sec- 
ondary tracks are therefore very likely to be due to aliasing with 
subtle effects on the count-rate or spacecraft orientation occurring 
on the Swift geo-orbital period. 

Figure [6] shows that ISGRI derived orbital solution is valid 
throughout the outburst and can be reliably extrapolated to times 
after the ISGRI observations ended at MJD 54576. The first BAT 
survey mode data available for the outburst, as shown in Fig. [T] 



is on MJD 54494; detector plane rate data is available before this 
time, although there does not seem to be any strong evidence for a 
pulsed signal before ~ MJD 54490. 

As is a tracer of the mass accretion rate, m, P can be used as 
an indicator of how m changes i.e. the rate at which the mass con- 
tained within the accretion disc is being exhausted. Consequently, 
we can define a simple timescale for the survival of the accretion 
disc, T, as: 



This gives t = 87 ±9 days based upon the ISGRI measurements 
shown in Table [2] Extrapolating from the time of peak of the out- 
burst this implies that the outburst should end at ~MJD 5461 1. 

4.5 Pulse Profile 

The shape of the pulse profiles has been measured over the 5- 
50 keV range with INTEGRAL. As seen in section |4T] the ISGRI 
timing results are highly accurate with strong power spectra leading 
to uncertainties on P of typically 200 /js. Phase diagrammes for the 
pulsar were made from the ISGRI 20-50 keV data from each revo- 
lution by folding on the measured period. An offset was then added 
to the phase diagram so that the brightest phase bin was centred 
at phase = 0. The period and epoch defined by ISGRI were then 
used to fold the JEM-X data in 5-25 keV band in order to make 
comparable phase diagrammes, which are shown in Figure|7] 

In general the shapes of the pulses seem to be more complex 
in shape below 20 keV than they are above 20 keV; below 20 keV, 
there are clear signs of a double peaked pulse profile. A secondary 
peak, approximately 50% the size of the primary peak, is especially 
noticeable in revolutions 654 and 655. In the same revolutions the 
pulse profile also shows signs of a secondary peak above 20 keV, 
although its relative size is much reduced. The pulse appears weak 
in revolutions 660-662; in revolution 660 the pulse can barely be 
detected at all. It is interesting to note that these three observations 
straddle the apastron of the system, when the separation between 
Be star and NS is largest and presumably the rate of accretion from 
any circumstellar material would be at a minimum. The suggestion 
is therefore that there is a connection between the orbital phase and 
the strength or shape of the pulsation. However, the observations 
of revolutions 670-672 are also taken at approximately the same 
orbital phases, and here a very clear sinusoidal pulse is visible. 
Unfortunately, determining any connection with the behaviour of 
GRO J 1750- 27 and the orbital phase would certainly require ob- 
serving an outburst over a larger number of orbital cycles and is not 
possible here. The variation in the pulse shapes from observation 
to observation show that it would not be appropriate to attempt to 
improve the signal to noise by making an average of many pulses. 

4.6 Pulsed Fraction 

The e volution of the pulsed fraction as defined in iBildsten et al.l 
( I1997I) , namely the ratio of the integrated peak-trough flux divided 
by the mean (pulsed-l-unpulsed) flux, is shown in Fig. [8] The pulse 
profiles in Fig.Qwere rebinned by a factor of two in order to calcu- 
late these values. The pulsed fraction in the 5-20 keV band varies 
between 5-15 %; above 20 keV the pulsed fraction is higher and 
has a larger variation, ~ 10^0 %. Both the 5-20 keV and 20- 
50 keV data show similar evolutions of the pulsed fraction with 
a slow rise in the first few observations, followed by a minimum 
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Figure 6. Dynamic periodogrammes measured from 50 ksec sections of Swift/BAI detector plane count-rate data for tlie period MJD 54480-54640 (2008 
January 15 - June 23) using the Bayesian analysis method outUned in section l4!4l Pulsation probability is represented with a rainbow colour scale with red 
representing high confidence of a pulsed signal. The central dashed line shows the orbital solution for the pulse variation measured from the ISGRI data 
(see section I4~l1 : the upper and lower dashed lines show the model offset by ± 3.4 ms. The ISGRI orbital model is in veiy good agreement with the BAT 
observations, no scaling or normalising of the ISGRI model has been performed. The offset lines agree with the secondary tracks caused by the aliasing of the 
pulsed signal with the Swift spacecraft orbital period. 



at revolution 660. In the 5-20 keV band the pulsed fraction seems 
to recover quickly to around 12 % in revolution 661 and remains 
approximately constant from then on; this behaviour is quite differ- 
ent above 20 keV, where the pulsed fraction increases slowly from 
revolution 661 towards the end of the INTEGRAL observations. 



5 DISCUSSION 

The pulsar spin-up and orbital parameters of the BeXRB system 
GRO J1750-27 have been measured to a high precision from hard 
X-ray observations with INTEGRAL and Swift (see TableO. These 
w ere found to be in good agreement with the measurements made 
bv lScottetal1 i ll997l) during the previous known outburst from this 
system in 1995. Assuming a cut-off power law model to describe 
the spectrum measured with INTEGRAL during the peak of the out- 
burst and extrapolating to 0.1-100 keV gives a total flux from the 
source of 6. 5x10"'' ergs cm"- s"'; assuming that during the peak 
of the outburst the source is radiating at the Eddington luminosity, 
l.SxlO^^ergs s"', puts an estimate on the distance to the source of 



around 15 kpc although this value is very uncertain due to the dif- 
ficulties measuring and then extrapolating the 5-50 keV spectrum 
and also assumes spherically symetric accretion onto the NS. 

The observed rapid spin-up of X-ray pulsars in outburst is pro- 
portional to the rate of angular momentum transfer from the donor 
star, which is in turn proportional to the rate of matter transfered 
to the NS through accretion, and also proportional to the X-ray lu- 
minosity of the source (e.g |Pringle&Reelll972h . In binary sys- 
tems with an accretion disk, the rate of spin-up is curtailed by 
the threading of the accretion disk with magnetic field lines from 
the NS, which causes a torque to be imparted to the NS and acts 
to reduce the rate of spin-up. This model was first proposed by 
iGhosh & Laiiibl j 19791) and alt hough other more detailed models 
have been proposed (e.g'Wang' 1987, 1995') they are largely con- 
sistent with the original. Ghosh & Lamb ( 1979) suggest that for a 
NS with radius 1 x 10* cm and given magnetic moment, yu, the 
value of P is dependent only on the quantity PL^^'', where L is the 
bolometric luminosity of the source. The changing luminosity can 
be estimated from the SwiftfBKT light curve (fig[T)> which has ex- 
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Rev 651: MJD 54507.687: Period 4.45450s: Orbital Phase 0.55 Rev 654: MJD 54516.572: Period 4,45270s: Orbilal Phase 0.85, Rev 655: MJD 54519.562: Period, 4,451 75s: Orbital Phase 0.95 








Rev 657: MJD 54525.545: Period^ 4.45515s; Qrbitol Phase 0.15^ Rev 658: MJD 54528.538: Period 4,4535 1 s; Orbilal Phose 0.25^ Rev 659: MJD 54531.530; Period^ 4,45535s: Orbital Phase 0.55^ 





Rev 661: MJD 54557.773: Period 4,45280s; Orbilal Phase 0.56, Rev 662: MJD 54540.509; Period, 4,4520Bs; Orbital Phase 0.65 



Rev 665: MJD 54545.752; Period^ 4 451B0s; Orbital Phase 0.76^ Rev 664: MJD 54547.759; Period^ 4,450585; Orbilal Phase 0.90^ Rev 667: MJD 54557.708; Period^ 4,4520Bs; Orbital Phase 0.23^ 





Rev 670: MJD 54564.432: Period 4.451935: Orbital Phase 0.46 Rev 671: MJD 54567.682: Period 4.45169s: Orbilal Phase 0.56 Rev 672: MJD 54570.413: Period 4.45130 : Orbital Phase 0.66 






Rev 675: MJD 54575.403; Period^ 4.450B1s; Orbital Phase 0.76 Rev 674: MJD 54576.595; Period^ 4.449965; Orbilal Phase 0. 





Figure 7. Relative pulse profile diagrammes for difi'erent INTEGRAL observations of the 2008 outburst of GRO J1750-27: ISGRI 20-50 keV (top, thickest); 
JEM-X 5-20 keV (bottom, thinnest). The data were folded on the periods listed in Table[T]and are indicated at the top of each plot along with the start date of 
the observation and the estimated orbital phase (0 = periastron) at that date. For clarity, the profiles are repeated for two phases. The numbers to the left of the 
profiles indicate the pulsed fraction and uncertainty. 
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cellent time coverage of the outburst and consists of average flux 
measurements taken daily. 

The BAT count rate was converted into erg cm^^s"' by 
comparing with the many Swift/BKI observations of the Crab 
Nebula; the average Crab Nebula count rate in 15-50 keV is 
0.22 counts cm^^s"'. In the same energy range, the Crab nebula 
was estimated to be 1.6 x 10~^ erg cm~^s~', assuming a photon 
spectrum of dN/dE = 10£"^°^ ("Zombeck' 2007). In order to find a 
bolometric luminosity, it is assumed that most of the energy output 
from GRO J1750-27 is in the range 0. 1-100 keV. The spectrum de- 
scribed in section[3]returns a 15-50 keV flux that is approximately 
a factor of 3 less than the total flux. Hence the BAT count rates can 
be converted to total energy output in ergs by scaling by the factor 
3 X (1.6 X 10"*'/0.22). This leads to a simple expression for the lu- 
minostiy in units of 10'^ ergs s"', L ^ 2.6F^p^d^, where Fbat is 
the BAT measured count rate in the 15-50 keV band and d is the 
distance to the source in kpc. 

The intrinsic spin-up model (Equation[5]l was then used to cal- 
culate P for the duration of each daily BAT measurement. The re- 
sults, calculated for v arious values of d, are presented in Fig. |9] 
against the models of lOhosh & LambI ( 1 19791) for a 1.4 M© NS 
with a radius of 10* cm and various values of ji. Several results 
are shown for the BAT data scaled for different distances. 

It can clearly be seen in Fig. |9] that from the peak of the 
outburst, when PLr'l'' is highest, that the BAT data show a slope 
that closely matches the gradient of the model, before turning over 
when the magnetic torque is expected to start to brake the pulsar. 
This requires confidence in the shape of the BAT light curve; a 
similar effect could be caused if there was an over estimation of the 
flux from GRO J1750-27 towards the end of the outburst. In sec- 
tion |4.4| the spin-up of the pulsar was found to end at MJD 54611, 
which is when the 15-50 keV flux from GRO J1750-27 is approx- 
imately equal to that from GX 3+1 (Fig.[T](. The separation of the 
two sources is 18', sli ghtly less than the FWHM of the BAT an- 
gular resolution of 22' jBarthelmv ;t al. 2005). Assuming a Gaus- 
sian point spread function, where the FWHM is equal to 235a, the 
maximum possible contamination of a measured flux caused by an 
equally bright source located 18' away is ~ 5%. If the measured flux 
from GRO J1750-27 was over-estimated by this amount it would 
have a negligible effect on the appearance of the data in Fig. [9] - 
equivalent to causing a value of log(PL'''') = 1 to be changed to 
~ 0.97. In addition, the flux from GX 3+1 shown in Fig. [T] re- 
mains reasonably constant throughout the outburst and no signif- 
icant correlation was found between the two light curves. Finally, 
the ability of the BAT batcelldetect software, which is used to pre- 
pare the transient monitor light curves, to succesfully resolve two 
sources is thought to be accurate where the separation is more than 
10' (Rrimm, private communication). Therefore any contamination 
of the GRO JI750-27 flux is likely to be considerably less than 5%. 

The value of p which seems to most closely resemble the 
model is ~ 3 x 10^" G cm', corresponding to a magnetic field at 
the surface of the NS of ~ 2 x lO'^ G. Fig. |9] suggests a distance to 
GRO J1750-27 of around 16 kpc. However, estimating the distance 
from this model is uncertain due to the difficulties in estimating the 
spectrum of GRO J 1750-27 and the crude scaling of the flux val- 
ues against the Crab Nebula, which is quite a different object spec- 
trally to GRO J1750-27; a variation in the bolometric correction 
of a factor of 2, for example, would lead to a change in the dis- 
tances by a factor of V2. Nevertheless, this makes it very likely that 
the object is located some distance beyond the Galactic centre (ap- 
prox 8.5 kpc from Earth), which is reinforced by the non-detection 
of any optical counterpart. A value of 30 kpc, which would put 



GRO 11750-27 outside the Galaxy, seems unlikely; additionally, a 
distance greater than 16 kpc would mak e GRO JI750-27 the fur- 
thest HMXB yet detected by INTEGRAL ( Bodaghee e t al.l2007h . A 
reasonable estimate on the range of acceptable distances measured 
here is therefore 12-22 kpc. HMXBs are relatively young systems 
that should be found near to their birth places in star formation 
re gions in t he sp iral arms of the Galaxy; the 4-arm spiral model 
of lRusseill ( I2OO3I) gives distances to the Perseus (~ 14 kpc) and 
Cygnus (~ 16 kpc) arms that are within the range mea sured here. 
Note t hat the 18 kpc distance lower limit estimated bv lScott et al.l 
( Il997l) was based on an assumed pulsed fraction of 30%; if the 
pulsed fraction during the 1995 outburst was ~ 15%, as measured 
in section l46l their distance would be reduced to > 13 kpc. Whilst 
the distance estimates are u ncertain, the change in P due to the ac- 
cretion torque predicted bv lGhosh & LambI jl979l) seems to agree 
with the data. 

The strong magnetic field estimated at~ 2xl0'^G 
opens the possibility to observe a Cyclotron Resonance Scatter- 
ing Feature (CRSF) in the X-ray spectrum of GRO JI750-27. 
To date CRSFs have been detected from 16 accreting pulsars 
jSchonherr et alj2007h and in V 0332+53 h ave also been o bserved 
Krevkenbohm et al. 2005 ; Mowla vi et al.l 
I hooeh . In the case of GRO JI750-27 the 



to evolve with time 



I2OO6I ; iTsygankov et al.l 



fundamental CRSF should be centred near to 22 keV; this would 
be challenging to observe with INTEGRAL as the peak would lie 
near to the boundaries of the ISGRI and JEM-X sensitivities and 
the peak flux from GRO J1750-27 is about a factor of 5 less than 
for V 0332+53, but is something that will be investigated in the fu- 
ture. Similarly, the monitoring capabilities of BAT could allow very 
detailed measurements of the evolution of CRSFs during a long 
outburst, if the sensitivity in narrow spectral bins was optimised. 

It is interesting that in both 5-20 and 20-50 keV energy bands 
the pulsed fraction maxima during the early part of the outburst 
appear near the periastron of the binary system, and that the minima 
are near apastron; however, the connection between the bands is 
broken in INTEGRAL revolutions 661-664 and in the following 
orbit no rise and fall of the pulsed fraction is visible within the 
error bars. Unfortuantely, even if there where some correlation of 
the behaviour of GRO JI750-27 with the orbital period it would 
be very difficult to confirm it with observations that sparcely cover 
only a few orbital phases. The behaviour of the pulsed fraction with 
the overall brightness of the source is shown in figure [10] Above 
20 keV there is a clear evolution of the pulsed fraction with flux, 
with opposite behaviour below and above 30 cps: at the beginning 
of the outburst (right of figure |I0| ) the pulsed fraction decreases 
proportionally as the source dims towards 30 cps; below 30 cps the 
pulsed fraction increases with decreasing flux. Note that the first 
point above 30 cps is the one taken during revolution 660, where 
the pulsation was at its weakest in 20-50 keV and not measurable 
in 5-10 keV (see figure[7]l. This also seems to coincide with a large 
drop in the average 20-50 keV flux between revolutions 660 and 
661. The same effect is not visible in the lower energy band which, 
except for a minimum again during revolution 660, shows a slow 
increase in the pulsed fraction with reducing count rate. 

The different behaviour above and below 20 keV rules out the 
attenuation of the pulsed flux by some intervening material that 
slowly accumulates and then dissipates; in this scenario the lower 
energy X-rays would show at least as strong an effect as at high 
energies. Work continues to investigate these effects with fine res- 
olution temporal-spectral analysis. 

On their own, the pulse period measurements made using the 
BAT non-imaged count-rate data would be difficult to interpret be- 
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Figure 8. INTEGRAL measurements of the pulsed fraction during the 
outburst. The dashed vertical lines show the times of periastron of the 
GRO J 1750-27 system. 
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Figure 9. Plot of P, taken from the spin-up model measured based on the 
INTEGRAL data (section l431 . against the quantity Pl}^"^ derived from the 
15-50 keV 5vvr/f/BAT count rates scaled for various distances (points). Each 
series of points is labelled with the distance, in kpc, used to estimate the lu- 
m inosity from the Swi/f/B AT count rates. The different lines show the model 
of iGhosh & LambI h979l) for a 1.4 Mq, \& cm NS for vaiious values of 
the neutron star magnetic moment, jU (indicated towards the bottom of each 
line). The units of P, P, L and /j are s yr s, 10^' ergs S-' and lO^" G cm^ 
respectively. The value of log|Q(PL-''^) for a 4.45 s pulsar emitting at the 
Eddington luminosity of ~ 2 X 10^^ erg s~' is approximately 1.2; a value of 
1.4 is approximately 3 times LecU- 
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Figure 10. Plots of the INTEGRAL measurements of the pulsed frac- 
tion against total pulsed+unpulsed flux during the 2008 outburst of 
GRO J1750-27. 



cause of the ambiguities in identifying the correct period, due to 
aliasing of the periodic signal, and the poor quality of the data 
caused by the background sources in the very large field of view. 
However, once these results are compared with the very high qual- 
ity pulsation measurements from ISGRI it is very easy to see that 
the time evolution of the pulsar can be measured over a long time 
scale with BAT. This suggests that outbursts from other BeXRBs 
and other reasonably bright pulsars with periods > few seconds 
could be equally well measured and shows the power of having 
two complimentary instruments in Syvift and INTEGRAL simulta- 
neously operational. If a proper treatment of the background in the 
detector plan count-rate data led to an increase in sensitivity then 
it may be feasible to follow the pulsations from numerous objects, 
for example the large numbers of BeXRBs found in the Magellanic 
clouds. 
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